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Abstract 

Sensor network localization problem is to determine the position of the sensor 
nodes in a network given pairwise distance measurements. Such problem can 
be formulated as a polynomial minimization via the least squares method. 
This paper presents a canonical duality theory for solving this challenging 
problem. It is shown that the nonconvex minimization problem can be re- 
formulated as a concave maximization dual problem over a convex set in a 
symmetrical matrix space, and hence can be solved efficiently by combining 
a general (linear or quadratic) perturbation technique with existing opti- 
mization techniques. Applications are illustrated by solving some relatively 
large-scale problems. Our results show that the general sensor network lo- 
calization problem is not NP-hard unless its canonical dual problem has no 
solution. 

Keywords: Sensor network localization; Canonical duality theory; 
Perturbation method Global optimization; NP-Hard problems 



1. Introduction 

Sensor network localization |3j,j6[, 18|, 26|, 34| is an important problem 
in communication and information theory, and has attracted an increasing 
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attention. The information collected through a sensor network can be inter- 
preted and relayed far more effectively if it is known where the information is 
coming from and where it needs to be sent. Therefore, it is often very useful 
to know the positions of the sensor nodes in a network. Wireless sensor net- 
work consists of a large number of wireless sensors located in a geographical 
area with the ability to communicate with their neighbors within a limited 
radio range. Sensors collect the local environmental information, such as 
temperature or humidity, and can communicate with each other. Wireless 
sensor network is applicable to a range of monitoring applications in civil and 
military scenarios, such as geographical monitoring, smart homes, industrial 
control and traffic monitoring. There is an urgent need to develop robust and 
efficient algorithms that can identify sensor positions in a network by using 
only the measurements of the mutual distances of the wireless sensors from 
their neighbors, which is called neighboring distance measurements. The ad- 
vance of wireless communication technology has made the sensor network a 
low-cost and highly efficient method for environmental observations. 

Sensor network localization can also be formulated as an optimization 
problem by least squares method. However, this problem is nonconvex with 
many local minimizers. To find global optimal solutions by traditional the- 
ories and local-search methods is fundamentally difficult. It turns out that 
the general sensor localization problem has been considered to be NP-hard 
[2I, |23[. Several approximation methods have been developed for solving this 



difficult optimization problem (see, for example, j2J|, j25|, j30|], [32| and refer- 



ences cited therein). The semi-definite programming (SDP) and second-order 
cone programming (SOCP) relaxations are two of the most popular methods 



studied recently j^, |25|, [30|, |32|. The basic idea of SDP relaxation is to 
think of the quadratic terms as new variables subject to a linear matrix in- 
equality. The SOCP relaxation is developed in a similar way. For both SDP 
and SOCP relaxation, computed sensor locations are not accurate when the 
solution of the localization problem is not unique. This is because many 
numerical schemes, such as primal-dual and interior point methods, for SDP 
or SOCP relaxation often return to the analytic centre of the solution set. 
These solutions are, in general, not global optimal solutions. 

Mathematically speaking, the localization problem in can be stated 
as follows [ij, jsf: Consider a sensor network in with m anchors and n 
sensors. An anchor is a node whose location G Mf^, where k = 1, ■ ■ ■ ,m, 
is known, and a sensor is a node whose location Xi G M"', where i = 1, ■ ■ ■ ,n, 
is yet to be determined. For a pair of sensors Xi and Xj, their Euclidean 
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distance is denoted as dij. Similarly, for a pair of sensor Xi and anchor a^, 
their Euclidean distance is denoted as Cik- In general, not all pairs of sensor 
/sensor and sensor/anchor are known, so the known pair- wise distances of 
sensor/sensor and sensor/anchor are denoted as G Ad and {i,k) G Ae, 
respectively. However, if we directly apply the general least squares method, 
the computation is very expensive and not practical for large problems [l7| . 

Canonical duality theory developed from nonconvex analysis and global 
optimization (see |8|, [isj) is an effective methodology, which has been used 



successfully for solving a large class of challenging problems in various dis 



ciplines. See, for example, [10|, ll2|, ll3|, |33|. This presents an effective per- 



turbation method based on the canonical duality theory to solve the general 
sensor network localization problem. Our main contribution is to show that 
this nonconvex optimization problem is not NP-hard unless its canonical dual 
problem has no solution. The rest of this paper is organized as follows. In the 
next section, we first rewrite the original problem as a new problem, where 
the decision variable is expressed in tensor (matrix) forms. In Section 3, the 
canonical duality theory is discussed in tensor space and a general analytical 
solution form is obtained by a complementary-dual principle. In section 4, 
the general sensor localization problem is first reformulated in vector space 
and then transformed as a concave maximization dual problem over a con- 
vex feasible space S^. Based on the triality theory, a quadratic perturbation 
method is proposed, which shows that the nonconvex sensor network opti- 
mization problem is not NP-hard unless its canonical dual problem has no 
solution in S^. Section 5 presents some concrete numerical experiments for 
sensor localization problems with two, 18, 20 and 200 sensors. The cases with 
noise are also considered. Results are compared with standard semi-definite 
programming method. Concluding remarks are given in the last section. 

The notations used in this paper are: M denotes the set of real numbers; 
A^denotes the transpose the matrix A. For a finite set S, \S\ denotes its 
cardinality and the bilinear form {u, u*) is simply the scalar product of two 
vectors or tensors. 



2. Problem Statement 

We use the least square method for the formulation of a new optimization 
problem. Consider a general sensor network localization problem, where the 
sensor locations are to be determined by solving the system of nonlinear 
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equations 



[i,]) e Ad, 

{i,k) G Ae. 



(1) 

(2) 



Here, the vectors a/^, A: = 1, • • • , m, are specified anchors, where 



I I 



\ a=l 



denotes the Euclidian distance between locations Xj and Xj G M.'^, i 



1, ■ ■ ■ , n, and 



= {(hj) e [n] X [n] 
Ae = {{i, k) G [n] X [m] 



|x,; — Xjll = (ijj, i < j, dij are given distances}, 
||xj — a^ll = Cik, Cik are given distances}, 



where [A^] = {1, ■ ■ ■ , N} for any integer N. 

For a small number of sensors, it might be possible to compute sensor lo- 
cations by solving equations ([I])-(l2l). However, solving this algebraic system 
can be very expensive computationally when the number of sensors is large. 

By the least squares method [2^, the sensor network localization problem 
(Vo) can be reformulated as a fourth-order polynomial optimization problem 
stated below. 



Wij[\\yii - Xji 



CP,): min | n(X) = J] 1 

{i,k)eAe J 



4) 



(3) 
(4) 



.Xn] = {Xf} G 



pdxn 



where X = [xi, X2, ■ ■ ■ 
Xi being a position in 
true sensor locations if and only if the optimal value is zero. This noncon 



is a matrix with each column 
Wij,qik > are given weights. Obviously, X are 



vex optimization problem appears extensively in mathematical physics [16 



computational biology [33(], numerical algebra [27| as well as finite element 
analysis of structural mechanics 28|]. Due to the nonconvexity, this prob- 



lem could have many local minimizers. It is fundamentally difficult, or even 
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impossible, to find global optimal solutions by traditional direct methods. 
In the following, we shall see that by using the canonical duality theory, 
this nonconvex minimization problem is equivalent to a concave maximiza- 
tion dual problem over a convex set, which can be solved efficiently by a 
proposed perturbation method. 

3. Canonical duality theory: A brief review 

The canonical duality theory is composed mainly of 1) a canonical trans- 
formation; 2) a complementary-dual principle; 3) a traility theory. This the- 
ory can be demonstrated by solving the following general nonconvex problem 
(the primal problem {V) in short) 

(V) : min (n(x) = ^(x. Ax) - (x, f ) + W-(Bx)| , (5) 

where Xa C R"^^" is a given feasible space, (x, x*) denotes the bilinear form 
between x and its dual variable x*, f G X* C R"^*^ is a given matrix, 
A : A'a — )■ X* is a given self-adjoint linear operator, B is a linear operator 
which assign each x G A'a to a (deformation gradient-like) variable in a 
linear space Wa, on which, W(w) : Wa — t- R is a well-defined differentiable 
nonconvex function. 

The canonical transformation is to choose a "geometrically admissible" 
nonlinear operator (see [sl) 

^ = A(x) : X, ^ S,, (6) 

which maps the convex set 3^^ into a convex set Sa, and a canonical function 
V : £a ^ such that the nonconvex functional W{w) can be recast in a 
canonical form 14^(Bx) = V{A{x.)). Thus, the primal problem {V) can be 
written in the following canonical form: 

(V) : min {n(x) = V(A(x)) - [/(x)} , (7) 

. r-i 

where f/(x) = (x, f) — ^{'x, Ax). By the definitions introduced in |8i], a non- 
linear operator A(x) : Xa Sa is said to be geometrically admissible if it 
can be used as a (deformation) measure such that the canonical transforma- 
tion 14^ (Bx) = V^(A(x)) satisfies certain necessary (geometrical and physical) 
conditions, such as the objectivity and isotropy. Let 

7^ = {Rg R™^"| = R-\ detR=l} 
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be a proper orthogonal rotation group in R*". 

Definition 1 (Objectivity and Isotropy |@] ). ^ subsetWa is said to be 
objective if 

Rw eWa Vw G Wa and VR G 7^. (8) 

A real-valued function W : Wa — M zs said to be objective if its domain is 
objective and 

W{Rw) = W{w) Vw G Wa and VR G 7^. (9) 

A subset Wa is said to be isotropic if 

wR G Wa Vw G Wa and VR G 7^. (10) 

j4 real-valued function W : Wa — )■ M «s said to be isotropic if its domain is 
isotropic and 

W{wR) = W{w) Vw G Wa and VR G 7^. (11) 

Geometrically speaking, the objectivity means that the function Vr(w) 
does not depend on rotation, but only on certain objective measure of its 
variable w. Therefore, the most simple objective function is the right Cauchy- 
Green deformation tensor C = w^w since 

C(Rw) = w^R^Rw = w^w = C(w) ^ VR G 7^. 

While the isotropy implies that the function W^(w) possesses a certain sym- 
metry. Clearly, the left Cauchy- Green deformation tensor ww^ is an isotropic 
measure due to the fact 

(wR^)(wR^)^ = ww^ ^0 VR G 7^. 

The concepts of objectivity and isotropy play important role in Semi-Definite 
Programming (SDP) and integer programming jo], [ll[. Particularly, if w is 
a vector, the objectivity is identical to isotropy. Furthermore, if the objec- 
tive function W(w) is considered as a kinetic energy and t/(x) is viewed 
as the potential energy, then the function 11 (x) is the true Lagrangian in 
mathematical physics. 

The objectivity in science is also refereed as frame invariance, which 
lays a foundation for mathematical physics and systems theory. In fact. 
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the canonical duality theory was originally developed from this concept |8|, 
which is the reason why this theory can be applied not only for modeling 
and analysis of complex systems, but also for solving a large class of noncon- 
vex/nonsmooth/discrete problems in both mathematical physics and global 
optimization (see review article fl3|). 

A differentiable function V{$,) is said to be a canonical function on its 
domain £a if the duality mapping <j = W{$,) from £a to its range S* is 
invertible. Let ^) denote the bilinear form on £a x £a- Thus, for the 
given canonical function V^(^), its Legendre conjugate V*{^) can be defined 
uniquely by the Legendre transformation 

y*(<j) = sta{(|;<j)-V(^) I ^ gU, (12) 

where the notation sta{g{^)\ ^ G Sa} stands for finding stationary point of 
g{^) on Sa- It is easy to prove that the following canonical duality relations 
hold on Sa- 

= VV(^) ^ ^ = Vr*(<?) ^ r(0 + V*i<;) = (13) 

By this one-to-one canonical duality, the nonconvex term 14^ (Dx) = V(A(x)) 
in the problem (V) can be replaced by (A(x); — V*{^) such that the non- 
convex function n(x) is reformulated as the so-called Gao and Strang total 
complementary function [8j: 

S(x,<j) = (A(x);<^)-V*(<^)-f/(x). (14) 

By using this total complementary function, the canonical dual function 
n°'(<j) can be obtained as 

n'^(<?) = sta{H(x,<j) I xG A-J 

= U\,)-V*{^), (15) 

where [/^(x) is defined by 

f/^(<?) = sta{(A(x); <j) - f/(x) | x G Xa}. (16) 

In many applications, the geometrically nonlinear operator A(x) is usually a 
tensor-valued quadratic function 

A(x) = I i (x, H«x) I : M''^" ^ M"><", (17) 
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where H^/ : Xa X* \/k, Z = 1, . . . , n is a symmetrical linear operator. In 
this case, the canonical dual variable ^ G C M"^" is a symmetrical tensor 
and the total complementary function H can be written in the following form 

H(x, = i(x, G Wx) - V*(<^) - (x, f), (18) 

where 

G(<;) = A + ^^HHfci. 

k,l 

For any given ^ E S*, the criticality condition Vx2(x, <j) = leads to the 
canonical equilibrium equation G(<^)x = f. Let 

Sa = {^e£:\ detG(.?)^0}. (19) 

Then on Sa, the solution to canonical equilibrium equation can be written as 
X = G^^(<;)f. Therefore, replacing the primal variable x by this generalized 
solution in S, the canonical dual function (fTSll can be explicitly formulated 
in the form of 

U'{,) = -^-{f,G-\,)f)-V*{,). (20) 

Theorem 1 (Complementary-Dual Principle [s*]). Problem (V^) is canon- 
ically dual to (V) in the sense that if ^ is a critical point o/n'^(<j), then 

x = G-i(^)f (21) 

is a critical point of n(x) and 

n(x) = s(x,^) = n^(^). (22) 

Conversely, if 5c is a solution to {V), it must be in the form of l[21\) for critical 
solution q o/n°'(<j). 

This theorem has extensive applications in nonconvex analysis and global 
optimization 10|. In finite deformation theory, this complementary-dual 
principle solved a 50-years open problem 2^- Note that the feasible set Sa 
is not convex, in order to identity the extremality property of the critical 
solutions, we need to introduce the following subsets of Sa- 

5+ = G Sal G{<i) y 0}, 5- = G Sal G{<^) ^ 0}. 
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Theorem 2 (Triality Theory). Suppose the (x, ^) is a critical point of 
S(x, <j). The critical solution ^ is a unique global minimizer of {V) if and 
only if ^ is a global maximizer o/n'^(<j) on , i.e. 

n(x) = minn(x) ^ maxW^^q) = n'^(^). (23) 

If ^ & , then ^ is a local maximizer ofIl'^{<i) on its neighborhood So C S~ 
if and only if 5c is a local minimizer of (V) on its neighborhood Xo C Xa, i.e. 

n(x) = maxn(x) ^ maxn°'(<j) = n'^(^). (24) 

If ^ E and dimXa = dim Sa, then ^ is a local minimizer o/n'^(<j) on 
its neighborhood So C S~ if and only if ^ is a local minimizer of {V) on its 
neighborhood Xo C Xa, i.e. 

n(x) = mmn(x) ^ minn'^((;) = n'^(^). (25) 

Proof. First, we need to prove the convexity of S^ and S~ . For any 
given <j^, G S^ , we have 

9G{^') yO, (1 - 9)G{<;^) y W9 e [0, 1]. 

Therefore, 

eG{^') + (1 - e)G{q^) = GiOq^ + (1 - d)<;') yO ^ee [0, 1]. 

This shows that S^ is convex. Similarly, we can prove the convexity of S~ . 

By the fact that the total complementary function S(x, <j) is a saddle 
function on Xa x S^, we have 

min n(x) = min max S(x, <;) = max min S(x, <j) = max n'^(<?). 

Therefore, by the complementary- dual principle, the critical point x G A'a is 
a global min of 11 (x) if and only if the critical point ^ is a global max of n'^(<?) 
on S^ ■ Since S(x, <j) is strictly convex in x and concave in q on Xa x Sa , its 
saddle point is unique. 

The statement (|24l) can be proved easily by the fact that 

maxmaxHfx, <j) = maxmaxSfx, <j) V(x, <j) ^ Xo x So C Xa x S~ . 
The statement fl25|) holds onl y u nder the condition dimAf^ = dim^a, 



which is a recent result proved in [15|. □ 
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Remark 1. The saddle min-max duality theorem / f^) was first proved by 



Gao and Strang in finite deformation theory lis], while the double min and 
double max duality statements were discovered in 1996. Based on this triality 
theory, the nonconvex minimization problem {V) is equivalent (only if ^ 
^) to a concave maximization dual problem 

iV") : max{n'^(<;) ^ G 5+} . (26) 

Although n'^(<j) contains an inverse matrix o/G(^), this canonical dual prob- 
lem can be solved easily by some well-developed nonlinear optimization tech- 
niques (see M [^/) as long as contains at least one critical point o/n'^(cr). 
Otherwise, a perturbation method will be discussed in the next section. 

To explain the theory, let us consider a very simple nonconvex optimiza- 
tion in M": 

min|n(x) = ia(^||xf -A)2-x^f Vx G M"| , (27) 

where a, A > are given parameters. The criticality condition VP(x) = 
leads to a nonlinear algebraic equation system in 

a(^||xf -A)x = f. (28) 

Clearly, to solve this nonlinear algebraic equation directly is difficult. Also 
traditional convex optimization theory can not be used to identify global min- 
imizer. However, by the canonical dual transformation, this problem can be 
solved completely and easily. To do so, we let ^ = A(x) = |||xp G M, 
which is an objective measure. Then, the nonconvex function W^x) = 

1^, /'I 11^112 



2 ^2 I 



x|p — A)^ can be written in canonical form V{^) = |a(^ — A)^. Its 
Legendre conjugate is given by V^*(<j) = ^a~^<;^ + A^, which is strictly con- 
vex. Thus, the total complementary function for this nonconvex optimization 
problem is 

S(x, = i||x||\ - ^a-'e - A^ - x^f . (29) 

For a fixed ^ G M, the criticality condition VxS(x) = leads to 

<:x - f = 0. (30) 
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Figure 1: Graphs of the primal function n(x) (sohd) and its canonical dual function n'*(^) 
(dashed). 

For each ^ 7^ 0, the equation f l5U]) gives x = f/<^ in vector form. Substituting 
this into the total complementary function S, the canonical dual function 
can be easily obtained as 



The critical point of this canonical function is obtained by solving the fol- 
lowing dual algebraic equation 



For any given parameters a, A and the vector f G M", this cubic algebraic 
equation has at most three real roots satisfying ?i > > <J2 > <^3, and each 
of these roots leads to a critical point of the nonconvex function P(x), i.e., 
Xi = f/^i, i = 1,2,3. By the fact that ?i G 5+ = G M K > 0}, ^2,3 e 
= {<;GM|<;<0}, then Theorem 2 tells us that xi is a global minimizer 
of n(x), X3 is a local maximizer of n(x), while X2 is a local minimizer if 
n = 1 (see Fig. 1). If we choose n = 1, a = 1, A = 2, and / = |, the 
primal function and canonical dual function are shown in Fig. [H where, 
Xi = 2.11491 is global minimizer of -P(x), ^1 = 0.236417 is global maximizer 
of n'^(^), and n(xi) = -1.02951 = n'^(^i) (see the two black dots). Also it 
is easy to verify that X2 is a local minimizer, while 0:3 is a local maximizer. 




{H(x,O|VxH(x,^) = 0} 
Ff 1 

-V-i«"V-A^, V^^O. 



(31) 




(32) 
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4. Application to Sensor Network Localization Problem 

Now let us tern our attention for solving the general sensor network op- 
timization problem (Pi). For convenience, we transfer variables from matrix 
to vectors, and let 

y = [xl ■ ■ ■ ■ ■ ■ xl^ - ■ ■ x^]^ G M"'^ : Locations of sensors (variables). 
Weights for the optimization problem {Vi), 

Q = [ill ■ ■ ■ Qlm ■ ■ - Qnl - ■ ■ qnm]'^ ^ I^"™ : 

Weights for the optimization problem (Pi), 

d d 



Sums of squares of anchors 

j2 .2 t2 .2 iT 



d = [dii ■ ■ ■ ■ ■ ■ d^^ ■ ■ ■ d^^ G M"" : Squares of distances between sensors. 



^ ~ [^11 ■ ■ ■ ^Im ■ ■ ■ ^nl ■ ■ ■ ^nm\ ^ 



Squares of distances between sensors and anchors. 
Then, Problem iVi) can be written in a vector form given below. 

{V): min i n(y) = i^^^. (y^D^^.y _ rf^ )2 

+ E [y^E^uy - 2Aly + Y.^at,f - e\ \ \ , 

(j,A:)GAe \ "=1 / J 

where Ajfc and a"^ are components of the anchors a^ obtained from the ex- 
pansion of ||xj — afclp in equation (jlj), Ei\^ G M"'^^"'' is a diagonal matrix 
defined by 



ik 




/ifc 




with /jfc G M'^^'^ being the identity matrix corresponding to sensor i and 
anchor fc, so that the (1,1) entry of /j^ coincides with the (i, fc) entry of -Ej^- 
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Similarly, D^j is an nd x nd matrix defined by 



















































. 


In 





















with liijjj, hjjji G IR'^^'^ being the identity matrices, so that the (1,1) entry 
of lii coincides with the entry of the matrix D^j. For Ijj, 1^, Iji, they 
are defined similarly. Let 



Qfe = Aifc(y) = y^^ifcy - 2A^y, 

where Aj^ and Ajj. are, respectively, geometrical operators from M"*^ into 



(33) 
(34) 



and 



£e = {e e R"^"| > 0}. 
By introducing quadratic functions Vij : ^ M and Vik : ^ M. such that 

1 



and 



^iK^^jfc) = stajeifcCTifc - 14fe(eifc) | Cife ^ Ve} 



(35) 



(36) 



+ (^e^, - j '^ik, {i, k) e A, (37) 



Problem (7^) can then be reformulated in the canonical form given below: 



(P): min^n(y)= J] V^.,(A,,(y)) + ^.^(A^.(y))l y e 



t,nd 
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Note that the function Vij{^ij) and Vik{eik) are both convex. Their duahty 
relations are given respectively as 

(^ij = VVij(^ij) = Wiji^ij - d^), (ij) G Ad, (38) 

and 

(Tik = ^V^kieik) = Qik [eik + ^Kkf - ' k) e Ae, (39) 

where <^jj and are dual variables. Let Sd be the range of the duality 
mapping = VVij(^jj), and let Se be the range of the duality mapping 
Cik = VVifc(ejfc)- Then, for any given <j G 5^ and cr G Se, the Legendre 
conjugate V*j and V^|, can be uniquely defined by 

V:^i,,,) = sta{e;Q, -y.,(6,) |6,GV4 (40) 

and 

^ifc(o"ifc) = sta{e^o-ife - 1/ife(eifc) | e^^ G Ve} (42) 

= l^k^ik + (^el - ^) ^ (43) 

where sta{} denotes finding the stationary point of the statement within {}. 
Clearly, (^, <j) and (e, cr) form a canonical duality pair (see j^). The following 
canonical duality relations hold on both Sd x Sd and Se x Se 

<^^J = vv^,(e.,) ^ 6, = vv^*(Q,)^e5Q, = v;,(6,) + ^*(^:')' 

cTifc = VVifc(ejfc) <^ Cjfc = VV;fc(criyt) ej^aik = Vik{eik) + V*f.{aik), 
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respectively. So the generalized complementary function ([s]) is given by 

{i,j)eAa {i,k)eAe 

{i,j)eAa (j,fc)GAe 

= E ((y"A,y)^.-(i<4 + <' 

{i,k)£Ae 
1 _ 

= ly^G(<j,a)y-FV)y-i(W-y(<joO 



(44) 



where s o t := [siti, ^2^2, ■ ■ ■ , Sn^n]"^ denotes the Hadamard product of any 
two vectors s, t G M", 



F(cr) 



'Y '^^k^ik ■ ■ ■ Y^ 2a^crifc ■ • • Y^ '^(^l^nk ■ ■ • y~^2afcr„fc 



fc=i 



k=l 



k=l 



k=l 



(45) 



G(<j, a) = 2(Diag (Fi(<j)) + Diag (F^ia)) + Gsi<;)), 



(46) 



with 



15 



Em 
k=l ^Ife 

Em 
k=l 

Em 
k=l ^nk 

Em 
k=l ^rik _ 



where the notation Diag {Fi) represents a diagonal matrix with Fu, i = 
1, ■ ■ ■ ,n being its diagonal elements. For a fixed ^ E Sd and cr e Se, the 
criticality condition VyH(y, ^, cr) = leads to the following canonical equi- 
librium equation: 



G(<j,o-)y-F(o-) = 0. 



(47) 



Substitute the solution of this equation into (jS]), the canonical dual function 
can be formulated as: 



n'^(<?,'T) 



_i(Q-i)^(o- o a) - d^<j + a^o- - e^a, 



where 



F((t) 



^nk ' ' / , -^^k^rik 



k=l 



and 



Q ^ 



=1 




k=l 


1 


1 


1 


Wll 


Win 




1 


1 


1 




Qlm 





k=l 



(48) 



T 
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Therefore, the canonical dual problem can be written in the form given below. 

{V") : sta{n'^(<j,o-)| qeSd^aeSe). 
By Theorem [T|, we have following results: 

Theorem 3. Problem (V^) is a canonical dual of the primal problem {V) in 
the sense that if &) is a critical point of (T"^), then 

y = G-\^,a)F{a) (49) 

is a critical point of {V) and 

U{y) = U\^,a). (50) 

Theorem 12] shows that the nonconvex primal problem {V) is equivalent 
to its canonical dual problem (V^) with zero duality gap, and the solution 
of {V) can be analytically expressed by pOl) in terms of the canonical dual 
variables. The global minimizer can be identified by the saddle min-max 
duality theorem fl23|) . In this case, the feasible space should be 

5+ = {(.?, (T)eSdxS,\ G(<?, 0-) y 0}. (51) 

Based on the triality theory, the nonconvex primal problem {V) is equivalent 
to the following canonical dual problem: 

(PlJ: max{n'^(.j,^T)| (<;,^)g5+}. (52) 

Theorem 4. // (^, a) G is a critical point of the canonical dual function 
n°'(^, a), then it is a unique solution ©/(P^^^), the vector y = G~"'^(^, ct)F(ct) 
is a unique global optimal solution to the nonconvex sensor network optimiza- 
tion problem {V), and 

n(y) = min n(y) = max {<; , a) = , &) . (53) 

Theorem |4] shows that the nonconvex sensor network problem {V) can 
be transformed to a concave maximization problem over a convex space S^, 
which can be solved easily if ^ 0. However, if n'^(<j, cr) has no critical 
point in S^, the primal problem {V) could be NP-hard, which is a conjunc- 
ture proposed in [9]. By the fact that is an open set, if the primal problem 
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has more than one global minimal solution, the canonical dual problem 11'^ (<j) 
may have no critical point in S^. In order to solve this problem, a quadratic 
perturbation method can be suggested as the following: 

min max |Sfc(y,<?,(T) = S(y,<;,(T) + ^Pfc||y-yfcf - (y,<5)]> (54) 

where > is a given parameter, y^ G M"'^ is a given vector (previous 
solution in iteration process), d G M"*^ is a given linear perturbation vector, 
S^^ is a relaxed canonical dual feasible space defined by 

5+ = {(<?, cT)eSdxS,\ G(.j, 0-) + /i J ^ 0}, (55) 

where < p^ < is given relaxation parameter. 

Remark 2. zs onr understanding that a NP-hard nonconvex optimization 
problem usually has multiple global minimizers. This is due to certain sym- 
metry of the modeling. The main idea of the linear perturbation is to destroy 
this symmetry such that the associated perturbed problem has a unique so- 
lution. This method is well-known in engineering mechanics but was first 
introduced to global optimization by the authors in 2008 for solving quadratic 
equations via least squares method [27]. Recently, this method has been used 
successfully to solve a NP-complete max-cut problem (see fs^ J). For complex 
systems, how to correctly chose the linear perturbation vector 5 is still an 
open problem. Therefore, some high-order perturbation methods proposed in 
2010 for solving certain NP-hard nonconvex optimization problems ll i J. By 
the fact that the perturbed total complementary function Sfc(y, <j, cr) is strictly 
convex m y G M"''^ and strictly concave in (<j, cr) on the closed convex set S^^, 
an effective canonical primal-dual algorithm can be developed for solving the 
saddle min-max problem [54\ ). 

5. Numerical Simulations 

In this section, we will first look at a simple case of a network with four 
anchors and two sensors. The locations of the anchors are known while 
the locations of the sensors are to be determined. The linear perturbation 
method will be used to show how the symmetry can be destroyed such that 
this small scale network localization problem can be solved nicely. We then 
move on to formulate more general sensor networks by randomly generated 
test problems. The so-called root mean square distance (RMSD) will be used 
to measure the accuracy of the estimated positions. 
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Figure 2: Sensor network with two sensors and four anchors. 



5.1. A four-anchor sensor network localization problem 

Consider the sensor network problem with two sensors and four anchors, 
as shown in Fig |2l Let xi = X2 = [x2,X2]"^ G denote the 

locations of the unknown sensors and let aj, i = 3,4, 5, 6, denote the locations 
of the four known anchors. The sensor network location problem is to solve 
the following system of nonlinear equations: 

(Vo) ||xi - X2II = di2, ||xi - asll = ei3, ||xi - a4|| = cm, 

||x2 - asll = 625, ||x2 - aell = 626, 

where 

dl2 = ei3 = ei4 = 625 = 626 = 2, 

ag = [-2, VS = 1.7321]^, = [-2, -Vs = -1.7321'^, 
as = [2, = 1.7321]^, ag = [2, -Vs = -1.7321]^. 

Let y = [x\ x\ x\ e 

Then, the fourth-order polynomial n(y) is given by 

n(y) = \{{x\-x\f + {xl-xlf-2'Y + \{{x\ + 2f + {xl-V-^f-2'f 

+ + 2)^ + {xi + - 2^f + \{{x\ - 2f + {x\ - - 2r 

+ ^((x^-2r + (x^ + y3f-2Y. 

To solve this symmetrical sensor network problem (Vq), we use a linear per- 
turbation method (see [27l |) 

{Vs) : min {n,(y) = n(y) - d^'y \ y G R"} , (56) 
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where 6 = [Si Sf > is a given perturbation vector. On the canonical 

dual feasible space Sa defined by 

<Sa = ?2, <^3, <^4, ^5) Ul + ?2 + ?3 7^ 0, ?i + <;4 + ?5 7^ 0}, (57) 

the canonical dual problem to the 5-perturbed problem (Vs) is 

{Vt): max{n^(0kG54, (58) 



where 



n^W = -^F,(<j)^G-i(<j)F5(<^)-(rfi2)\i + (ara3-(rfi3)')?2 
+ (a4 a4 - ((ii4)^)g3 + (a^as - (^25)^)^4 

I / T /J \2\ ^2 ^2 ^2 ^2 ^2 

+ (a6a6 - (c/26) )^5 - - - - 2*^4 - 2^5 > 



Fs{<;) = [6l + 2al<;2 + '2alc;s,Sl + 2a% + '2a%, (59) 
+ 2a^^4 + 24<^5, + 2a^?4 + 2a2<;5]^, (60) 



G(0 



2(^1 + ^2 + ^3) -2c;i 

2(^1 + ^2 + ^3) -2^1 

-2qi 2(^1 + ^4 + ^5) 

-2^1 2(^ + ^4 + ^5) 



.(61) 



Set 5 = [0.005, 0.005, 0.005, 0.005]^. Then, the canonical dual problem has a 
unique solution 

^ = [-0.0000, 0.0005, 0.0020, -0.0020, -0.0005]^. 

By Theorem m it follows that 



r-i -2 -1 -2ir 

[X]^, X]^, X2, X2J 



[-0.9994, 0.0002, 1.0006, 0.0002]^ 



is a global minimizer of n5(y). 
It is easy to verify that 



Us{y) = -4.1667 x 10"^ = nf(^). 
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By the fact that 

n(y) = 4.1667 X 10"^ 

the (5-perturbed solutions y can be considered as the global minimizer to the 
original problem [Vq) and we have 

||xi - X2II = 2.0000, ||xi - agll = 2.0001, ||xi - a4|| = 2.0005, 
||x2 - asll = 1.9995, ||x2 - agH = 1.9999. 



5.2. 18 sensors network localization problem with four anchors 

We now consider sensor network localization problem with 18 sensors. In 
this case, we have Problem (Pi) with d = 2. Define y = [x\, xf, - ■ ■ , x^, x^]^ G 
M^", and let Wij = qik = 1 in Problem (Pi). Here, we do not consider noise. 

The 18 sensors {xj = [Xj^,a;f] : i = 1, ■ ■ ■ ,18} are randomly generated 
in the unit square [-0.5, 0.5] x [-0.5, 0.5]. The four anchors (ai, a2, as, a4) 
are placed at the positions (±0.45, ±0.45). The distances d = {dij}, i = 
1, ■ ■ ■ , 18; j = 1, ■ ■ ■ , 18, and e = {e^fc}, z = 1, ■ ■ ■ , 18; = 1, ■ ■ ■ , 4, are 
computed as follows: 

7 II * * 1 1 1 1 * 11 

Clij W'^i 1 1 ' W'^i ^k\\ 

We now assume that the locations of the 18 sensors are unknown. They 
are to be determined by the approach proposed in the paper. The sequen- 
tial quadratic programming approximation with active set strategy in the 
optimization toolbox within the Matlab environment is used to solve the 
canonical dual problem. 

By Theorem m we obtain y = [xi, ■ ■ • ,Xi8]-^ with x^ = [x.} , x.'^]'^ , i = 
1, • • ■ ,18, which is a global minimizer of n(y). 

Furthermore, we have 

n(y) = 1.30 X 10"^ ~ 3.03 x 10"^ = n'^(^, a). 

This problem is also solved by the standard semi-definite programming (SDP) 
method. The RMSD obtained using the canonical dual method is 4.61 x 10~^, 
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Figure 3: Sensor network with 18 sensors by the canonical dual method. 

while the RMSD obtained using the standard SDP method is 4.45 x 10~^, 
where RMSD is the Root Mean Square Distance defined by 

1 " 

RMSD = (-V||x,-x,f)5, 

n ^-^ 

1=1 

which is to measure the accuracy of the computed locations. 
The computed results by the canonical dual method and the standard SDP 
method are plotted in Fig. |3] and Fig. |U respectively. The true sensor 
locations (denoted by circles) and the computed locations (denoted by stars) 
are connected by solid lines. O ur p rogram is implemented in the MATLAB 
environment, where SEDUMI |2l| is used as the SDP solver. From the 
results obtained, we see that, when there is no noise and the sensor size is not 
too large, both the canonical dual method and SDP method are very effective 
method for finding sensor locations. In particular, for the canonical dual 
method, all the stars are exactly located inside circles. The computational 
time for canonical dual method and SDP method are 0.61 seconds and 21.53 
seconds, respectively. 
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Figure 4: Sensor network with 18 sensors by the standard SDP method. 

5.3. A 20-sensor-network localization problem with distance errors 

A network of 20 uniform randomly distributed unknown points is gener- 
ated in the square area [0, 1] x [0, 1]. We assume 

If II — Xj\\ < radio range, a distance (with noise) is given between 
Xi and 

If ||xj — Xj\\ > radio range, no distance is given between Xi and Xj. 

4 anchors are located in [0,0], [0,1], [1,0] and [1,1]. The distances between 
the nodes are calculated. If the distance between two nodes is within the 
specified radio range of 0.4, the distance is included in the edge set for solving 
the problem after adding a random error to it in the following manner: 

dij = + A^(0,0.001)| 

where dij is the actual distance between the 2 nodes, and A^(0, 0.001) is a 
random variable. 

The computed results obtained by the canonical dual method and the 
standard SDP method |19j] are plotted in Fig. [5] and Fig. [6|, respectively. 
The true sensor locations (denoted by circles) and the computed locations 
(denoted by stars) are connected by solid lines. The computational time for 
canonical dual method and SDP method are 0.65 seconds and 27.71 seconds. 
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Figure 5: Sensor network with 20 sensors solved by the canonical dual method. 
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Figure 6: Sensor network with 20 sensors solved by the standard SDP method. 

5.4- A 200-sensor-network localization problem with distance errors 

A network of 200 uniform randomly distributed unknown points is gener- 
ated in the square area [0, 1] x [0, 1]. 4 anchors are located in [0,0], [0,1], [1,0] 
and [1,1]. For all sensors, the radio range = 0.3. The distance, including a 
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Figure 7: Sensor network with 200 sensors solved by the canonical dual method. 
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Figure 8: Sensor network with 200 sensors solved by the standard SDP method. 



random error, is generated in the following manner: 

dij = + A^(0,0.001)| 

The computed results obtained by the canonical dual method and the stan- 
dard SDP method are plotted in Fig. [7] and Fig. [8], respectively. Careful 
examination of the results obtained for the cases involving 20 sensors and 200 
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sensors, we observe that when noise is taken into consideration, the canonical 
dual method gives rise to much better solutions. In particular, if the level 
of noise or the sensor size is large, the standard SDP is usually having diffi- 
culty to finding the exact sensor positions. See Figure [5] and Figure [3 The 
computational time for canonical dual method and SDP method are 127.10 
seconds and 1088.70 seconds, respectively. 

6. Conclusion Remarks 

We have presented a solid application of the canonical duality theory 
for solving challenging sensor network localization problems. By using the 
complementary-dual principle, a general form of tensor-valued analytical so- 
lution is obtained in terms of the canonical dual variables. Our results show 
that a nonconvex minimization problem could be NP-hard if its canonical 
dual problem has no solution in S^. Based on a quadratic perturbation 
method, an effective canonical saddle min-max approach is proposed for solv- 
ing this challenging problem. Applications are illustrated by detailed analysis 
of a small size problem as well as some relatively large scale sensor network 
localization problems. Finally, the performance of the method proposed on 
problems with noisy distance data needs further investigation. 
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